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Molecules of bisthiolterthiophene have been ad- 
sorbed on the two facing gold electrodes of a me- 
chanically controllable break junction in order to form 
metal-molecule(s)-metal junctions. Current- voltage (I- 
V) characteristics have been recorded at room temper- 
ature. Zero bias conductances were measured in the 
10-100 nS range and different kinds of non-linear I- 
V curves with step-like features were reproducibly ob- 
tained. Switching between different kinds of I-V curves 
could be induced by varying the distance between the 
two metallic electrodes. The experimental results are 
discussed within the framework of tunneling transport 
models explicitly taking into account the discrete na- 
ture of the electronic spectrum of the molecule. 



I. INTRODUCTION 

Molecular electronics, understood as 'making 
an information processing device with a single 
molecule' requires synthesizing molecules with 
electronics functionnalities and connecting them 
together and to external electrodes. In the sev- 
enties, A. Aviram and M. A. Ratner and 
F. L. Carter j|] proposed challenging design of 
molecules analogous to diodes or triodes. Since 
then, molecules of these types have been synthe- 
sized, and some connection techniques have been 
developed. State of the art is the connection of a 
few molecules, and even of a single one, to con- 
ducting electrodes. This was first achieved using 
a Scanning Tunneling Microscope by positioning a 
metallic tip above molecules deposited on a con- 



ducting substrat£ 
ious molecules 
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Electronic properties of var- 
, and eventually single Cgo 
ones [ |l7|Jl^ have been investigated in this way. 
The conclusion reached from these experiments 
is that a molecule can have a non-zero conduc- 
tance determined by its molecular orbital struc- 
ture ]l8|-p0t. Using a STM to contact a molecule 
however suffers from intrinsic limitations, such as 
the asymmetry of the junctions and, at least at 
room temperature, the lack of the mechanical sta- 
bility necessary to maintain a stable chemical bond 
between the molecule and the tip. These limi- 
tations, together with the improvement of litho- 
graphic techniques, have prompted the emergence 
of complementary techniques. In the last three 
years, two alternative techniques have been pro- 
posed for contacting a molecule to metallic elec- 
trodes. The first one consists in fabricating a se- 
ries of metallic electrodes, to which the molecule is 
contacted, at the surface of a substrate. This tech- 
nique, so far restricted to molecules longer than 
at least 5 nm ||2l|], has been used, in particular, 
to investigate carbon nanotubes [^,^. The sec- 
ond one is to use a mechanically controllable break 
(MCB) junction |p^ , p5| . It consists in breaking 
a small metallic wire, introducing molecules with 
end-groups reactive to this metal into the gap, 
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and adjusting the gap between the two facing elec- 
trodes to a distance comparable to the length of 
the molecule in order to contact it. This technique 
combines the advantages of the previous ones: it 
allows to connect a short molecule between similar 
electrodes while maintaining a very high stability 
(down to 0.2 pm/hour for nanofabricated MCB 
junctions) ||2^. As demonstrated recently by Reed 
et al |2j], it opens up the possibility to measure 
electrical transport through a very few molecules 
chemically bound to two electrodes, and possibly 
through a single one -the actual number being how- 
ever very difficult to know . 

In the present paper, we describe the use of gold 
MCB junctions to investigate the electronic trans- 
port properties through bisthiolterthiophene 2.5"- 
bis(acetylthio)-5.2'5'.2" 

-terthienyl (T3) molecules. We show (i) that dif- 
ferent types of reproducible I-V characteristics can 
be obtained and (ii) that varying the distance be- 
tween the two electrodes induces a switch between 
different I-V curves. Finally, we discuss the in- 
terpretation of the experimental results within the 
framework of two differents models which both ex- 
plicitly include the discrete nature of the electronic 
levels of the molecule: (i) a coherent model which 
treats the molecule as a scattering impurity be- 
tween two metallic wires, and (ii) a sequential tun- 
neling model, in which the molecule is assumed 
to be weakly coupled through tunnel junctions to 
each metallic electrode. 



II. EXPERIMENTAL TECHNIQUES AND 
RESULTS 
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FIG. 1. Ideal Sample. A conjugated molecule is 
chemisorbed onto the gold electrodes via the thiolate 
terminal groups. 
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FIG. 2. Representation of the T3 molecule. Acetyl 
protecting groups are visible (circled) at each end. 
They were removed prior to the assembly process. 

Suspended metallic microbridges were fabri- 
cated as described in Ref. [Q. First, an insu- 
lating layer of polyimide PI 2610 from Dupont 
dc Nemours was spun on a polished phosphor- 
bronze substrate. Using standard e-beam lithog- 
raphy techniques, a metallic nanostructure with 
the appropriate shape was then deposited on a 
4 mm X 20 mm chip. In our case, the metallic layer 
consisted in a 100 nm thick gold layer, with an un- 
derneath 0.2 nm thick titanium adhesion layer and 
a 5 nm thick aluminum protection layer on top. Fi- 
nally, isotropic reactive ion etching of the PI layer 
produced a metallic bridge suspended between an- 
choring pads (Fig.^. The suspended length was 
of the order of 3 /im, and the central constriction 
was less than 100 nm wide. 



A. Sample fabrication 

The aim of the experiment is sketched on Fig.|l|. 
The conjugated molecule T3, to be connected to 
the electrodes of the MCB junction, was syn- 
thesized from terthiophene. A thiolate function 
was substituted at both ends of the molecule for 
its ability to strongly react with gold surfaces 
|p7| . The thiolate functions were protected by 
thioester formation with acetic anhydride in order 
to avoid successive oxidative oligomerization that 
would generate polydisulfides |^ in solution; the 
protecting acetate groups were removed just be- 
fore immersion of the gold electrodes into the T3 
solution (Fi^h. 




FIG. 3. Scanning electron microscope picture of a 
suspended junction before breaking. 

The chip was then mounted in a three point 
bending configuration. The chip, resting on two 
countersupports, was bent by pushing on its cen- 
ter with a driving rod, actuated through a coarse 
adjustment screw, until the bridge breaks, as 
indicated by infinite electrical resistance. The 
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molecules were then self-assembled onto the freshly 
broken gold electrodes by immersion of the sus- 
pended junction in a droplet of a 5 x 10^"* M so- 
lution of T3 in trichloro-l,l,l-ethane. The droplet 
was supplied by a Hamilton syringe mounted 0.5 
mm above the junction. In order to improve the 
self-assembly process of the molecules, the acetyl 
protecting group of the dithiol was hydrolized in 
situ by adding 0.1 % of dimethylaminoethanol 1 
min before the experiment. The gold electrodes 
were kept in solution for 1 min, then the solvent 
was evaporated. After that. Argon was continously 
flushed through the shielded sealed box containing 
the set-up. These conditions were chosen so as 
to hinder the formation of di- or polydi-sulfides. 
In a final step, the bridge gap was reduced using 
a piezoelectric fine adjustment of the driving rod 
until a non zero conduction was detected (Fig. 0). 





FIG. 4. Actual experiment: the top and bottom 
gold electrodes are first separated by breaking the junc- 
tion (left), T3 molecules are adsorbed onto them (mid- 
dle) and the electrodes are brought closer to allow I-V 
measurements (right) on a single or a small number of 
molecules. 



B. Experimental setup 

All electrical connections were filtered with low- 
pass RC filters. The junction was voltage biased, 
and the current was measured using an I/V con- 
verter. Typically, 512 points were collected for 
each I-V curve by sweeping the voltage from -2 
V to 2 V in 0.2 to 20 seconds. 



C. Control experiments and Calibration 

As a control experiment, we first recorded the 
I-V characteristics of a metal-air-metal junction. 
I-V characteristics are linear, and the variations 
of the conductance with the piezoelectric actu- 
ator elongation are shown in Fig.||. Assuming 
elastic deformation of the substrate and a barrier 
height of 1 eV for a gold-air-gold tunnel junction 



[ P9[ , the observed exponential dependence |^0[ pro- 
vides a calibration of the displacement ratio r be- 
tween the piezo elongation and the inter-electrode 
spacing. The obtained value r = 3.3 x 10~^ is 
consistent with the geometrical estimate ||2^ , ^ 
r = 6tu/L'^ = 3.75 x 10"^ where t = 0.3 mm is 
the substrate thickness, u = 3 /im is the distance 
between the anchoring points, and L — 1.2 cm is 
the distance between the two countersupports. 
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FIG. 5. Conductance of a metal-air-metal junction 
as a function of the variation of the inter-electrode 
spacing (dots). The origin of the horizontal axis is 
arbitrary. The dashed line correspond to the WKB 
exponential variation assuming a barrier height of 1 
eV(see text). 

I-V characteristics of a metal-air-metal junction 
after 1 min immersion in the pure solvent were also 
measured. I-V curves wer featureless and exhibited 
a linear behaviour at low bias. 



D. Measurements of Au-T3-Au junctions 

During a typical experiment on Au-T3-Au junc- 
tions, stability periods with a duration 1-20 min 
alternate with instability periods generally lasting 
a few minutes. This behavior was always observed 
on the four samples that have been measured. A 
series of I-V characteristics, recorded subsequently, 
are shown in Fig.^. Although different I-V charac- 
teristics could be observed, the reproducible ones 
are of one of the two types of I-Vs shown in Fig.0. 
Asymmetric I-Vs of type (a) are more often ob- 
served and more stable than symmetric ones of 
type (b). 
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FIG. 6. Consecutive single sweep I-V curves 
recorded at room temperature for a gold-T3-gold junc- 
tion. Curves are shifted vertically for clarity. 
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FIG. 7. Typical (a) asymmetric (solid line) and (b) 
symmetric (dashed line) I-V curves recorded at room 
temperature for gold-T3-gold junctions. Both curves 
were obtained by averaging over 5 voltage sweeps. 

Fig.^ shows typical asymmetric curves. 
The measured zero bias conductance of type (a) 
junctions is of the order of 10 nS. The asymmetric 



I-V characteristic is non-linear with step-like fea- 
tures, and the current increases linearly at large 
voltage. 

The measured zero bias conductance of type (b) 
junctions is larger, about 80 nS. The symmetric TV 
characteristic is also non-linear with smaller step- 
like features. At > 1 V, the current rises faster 
than linearly with V. 

Fig.|| shows a mechanically induced transition 
between different TV curves. A series of repro- 
ducible type (a') TVs was first recorded. The 
inter-electrode spacing was then reduced by 0.04 
nm approximately, and a series of type (a") TVs 
were then recorded. It should be noted that upon 
reduction of the gap size, the zero bias conductance 
first decreased from 13 nS to 6 nS. 
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FIG. 8. TV curves recorded (a') before (solid line) 
and (a") after (dashed line) reduction of the in- 
ter-electrode spacing by ca 0.4nm. Both I-V curves 
were obtained by averaging over 5 voltage sweeps. 

III. DISCUSSION 

Before discussing transport models, we first 
present the electronic properties of the isolated T3 
molecule and their modification upon adsorption 
onto the electrodes. 



A. Electronic properties of a T3 molecule 

Fig. ^ (triangles) shows the electronic spectrum 
of an isolated deprotected T3 molecule calculated 
using a standard extended Hiickel technique, start- 
ing from a forced planar geometry optimized at 
the AMI level Although properties of the 

T3 molecule adsorbed onto gold are not known. 
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except for the adsorption of T3 in a matrix of 
dodecanethiol |Q, other rigid rods a,w-dithiols 
have been shown to form assembhes onto gold, in 
which one thiol group binds to the surface while 
the other thiol group projects upward at the outer 
surface of the self-assembled monolayer Fur- 
thermore, organised self-assembled monolayers of 
bis-(2,2':5',2"-terthien-5-yl)disulfide have been de- 
scribed From these studies, it seems plausi- 
ble that the T3 molecule, once deprotected, reacts 
onto the gold electrodes with (i) the formation of a 
Au-S bond and (ii) the terthiophene moiety point- 
ing upward rather than laying on the surface of 
the electrode. One key question concerns the posi- 
tion of the Fermi level of the electrodes relatively 
to the molecular electronic levels. The latter are 
expected to shift so that the Fermi level of gold 
falls in the HOMO-LUMO gap, its position being 
determined by the amount of charge displaced on 
the molecule. It has been shown that the formation 
of the Au-S bond involves a formal negative charge 
transfer , from the metal to the molecule. The 
values of (5^ ranges from 0.2 to 0.6 electron for thi- 
olate adsorption on gold Jl5| , ^^ |^ . Using the 
method described in reference P^ , we found a po- 
sition of the Fermi Level of —10.5 ± 0.06 eV in the 
energy scale used for the molecule. The HOMO 
is thus closer to the Fermi level than the LUMO, 
with Ef - Ehomo = 0.2 ±0.06 eV |38). How- 
ever, given the approximations made in the calcu- 
lation and the uncertainty about the exact amount 
of charge transferred, it seems reasonable to keep 
this difference as an adjustable parameter within 
the [0 — 0.7eV] range, i.e. Fermi Level between 
HOMO and midgap. 




FIG. 9. Transmission function of the T3 molecule 
calculated by the ESQC technique. The triangles 
represent the energy level positions of the isolated 
molecule ( filled occupied MO's and open unoccupied 
MO's of the deprotected T3 diradical). HOMO and 
LUMO refer to the gap of T3 adsorbed on the two 
electrodes. The energy scale reference is arbitrary. 



B. Linear regime 

At low bias {V < O.IV), the I-V characteristics 
are linear. The measured zero bias conductance 
are respectively 80nS, 6nS and 13nS for the curves 
in Figjl and Fig.|. 

We first calculated the transport properties of 
the gold-molecule-gold junction by the Electron 
Scattering Quantum Chemistry ( ESQC) tech- 
nique This technique has proven to 
quantitatively account for the low bias conduc- 
tance of single adsorbed molecules as probed by 
STM jl^]. It treats the molecule as a defect 
which breaks the translational invariance of the 
metal, and therefore scatters incident electrons. 
The ESQC technique ignores both the electron- 
electron and electron-phonon interactions and ne- 
glects charging effects. It assumes that the scatter- 
ing is elastic because for molecule of small length, 
the tunneling time is shorter than the intramolecu- 
lar relaxation times. An extended Hiickel model is 
used to build up the matrix representation of the 
multichannel scattering Hamiltonian taking into 
account the complete chemical description of the 
electrodes and of the molecule. The calculated 
multichannel transmission coefficient T{E) of an 
electron at a given energy E is shown in Fig. ^. 
The linear conductance G of the metal-molecule- 
metal junction is then determined using the Lan- 
dauer formula p| , 

2e2 

G = ^T{Ef) 

where Ef is the Fermi level of the electrodes. 

The prediction for G thus depends on the es- 
timated position of the electronic spectrum rela- 
tively to the Fermi energy Ep, on the exact con- 
formation of the molecule in the junction and of 
the coupling of the molecule to the electrodes. The 
strength of coupling is determined by the length of 
the Au-S bonds. The S atoms were assumed to 
be adsorbed in a hollow site of the gold surface. 
A bond length of 1.905 A was used in the present 
calculation |3^. This is the shortest distance, we 
found in the literature. It thus provides an up- 
per bound for the coupling strength. Assuming a 
symmetric coupling at both ends of the molecule. 
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the calculated conductances for Ep — Ehomo = 
0.6(midgap), 0.2 and eV are 87, 585, 2306 nS, re- 
spectively. Although the order of magnitude of 
these values is comparable with the measured one 
G ~80 nS for type (b) junctions the discrep- 
ancy indicates that the coupling of the molecule to 
the electrodes is smaller than estimated. This find- 
ing can be compared with the results published re- 
cently by Emberly and Kirczenow . In that pa- 
per, it is shown that the experimental conductance 
of gold-benzene di -thiol-gold junctions, reported 
in ref. can be accounted for by a transport 
model based on the Landauer formalism provided 
that artificially large bond lengths (i.e. small cou- 
plings) are used. In our case the discrepancy be- 
tween theory and experience seems less important 
but could also be reduced using longer Au-S bonds: 
for example Go=87nS at midgap for Au-S=1.905A 
reduces to Go=20nS at midgap for Au-S=3A. 

C. Non linear regime 

We now discuss the step-like features found in 
the I-V characteristics at large bias voltage and 
their possible interpretations. 

At first, we can rule out the simplified Coulomb 
blockade model put forward by Reed et al to ac- 
count for their experimental observations in a sim- 
ilar experiment with benzene dithiolate [ p4| . In- 
deed, this model, which considers the molecule as 
a usual metallic island forming a small capacitor 
with each electrode, predicts a series of steps which 
are not observed. 

Consequently, we attribute the steps to the dis- 
creteness of the molecular levels [Q . We present 
and discuss below two models which both involve 
the discrete electronic levels of the molecule explic- 
itly. The first one is a coherent tunneling model 
derived from the ESQC technique, whereas the 
second is a sequential tunneling model. Basically, 
they differ by the fact that the electron resides (se- 
quential tunneling) or not (coherent tunneling) on 
the molecule during the transfer. 

Given the symmetry of the molecule, symmet- 
ric I-V characteristics are expected in the case of 
symmetric couplings. Thus type (b) characteris- 
tics only will be analyzed below. 

1. Coherent tunneling model 

The coherent tunneling model is based on an 
extension of the ESQC technique outlined above. 
It supposes that the coupling to the electrodes is 
strong, in other terms that the tmmcling time of 



an electron through the molecule is much smaller 
than the intramolecular vibronic relaxation time: 
i.e. the molecule is not charged by the tunneling 
process, the electron having no time to be com- 
pletely relocalized in the molecule after the initial 
tunneling step at the metal -molecule electronic 
contact. The current is calculated using psf : 

I{V) = 4 r T (E, V) (/ {E -^i,)-f{E- /i2)) dE 

(1) 

in which the effect of the bias voltage V on the 
T(E) spectrum is included; /ii.2 is the chemical 
potential of the electrode 1,2, and /(e) denotes 
the Fermi function at the temperature of the ex- 
periment. In the case of a symmetric junction, 
the voltage at the molecule lies halfway between 
the electrode voltages. In order to calculate the 
current, we make the crude approximation that 
T{E,V) ~ T{E - eV/2,0) where T{E,0) is the 
spectrum of Fig. ^. Equation |l| then predicts the 
I-V characteristic shown in Fig. |lo| . 
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FIG. 10. Experimental symmetric I-V curve (dia- 
monds). I-V curve (solid line) calculated with the co- 
herent model. The calculation (see text) has been done 
using: Ef — EHomo=0.4eV. 

For both polarities, the first resonances occur 
through the HOMO and the HOMO-1 levels. This 
results from the assumption that the equilibrium 
Fermi level is closer to the HOMO than to the 
LUMO ||l^. In this model, the shape and height 
of the steps is only due to the broadening of the 
molecular levels by the coupling with the elec- 
trodes. In particular, the height of a step in the 
TV curve is directly proportional to the area of 
the corresponding peak in the T(E,V) spectrum. 



6 



The experimental curve is well reproduced for the 
[— 0.5y, 0.5V] range. Outside of this range, the 
calculated current is higher than the experimen- 
tal one. It should be noted that the only ad- 
justable parameter used in the calculation is the 
position of the Fermi level. Here we have taken 
Ep — Ehomo — 0.4 eV. In this calculation, the 
magnitude of the current depends on the following 
points: i) the overlap between the orbitals of the 
outer Au atoms and of the sulfur atoms; ii) the lo- 
calisation of the molecular orbital and its symme- 
try, iii) perfect symmetry between the two contacts 
p6t . The applied voltage is expected , especially at 
high bias, to affect points ii) and iii), both of which 
should tend to decrease the calculated current and 
thus reduce the discrepancy. 



2. Sequential tunneling model 

In this model, the molecule is treated as a 
quantum dot with discrete energy levels weakly 
coupled to both electrodes through tunnel junc- 
tions. It considers the tunneling of an electron 
through the metal-molecule-metal junction as a 
sequential process |Q, the molecule being suc- 
cessively charged and discharged. At first sight, 
the existence of tunnel barriers at both ends of 
the molecule can be questionned since the sulfur 
species chemically adsorbed on the electrodes do 
contribute to the HOMO and LUMO, both or- 
bitals being delocalised over the entire molecule. 
However, the existence of barriers to the injection 
of carriers at metal-oligothiophene junction have 
been reported in the literature p8[ |. Furthermore, 
the possible existence of a barrier to injection at a 
gold -thioacetylbiphenyl interface has also been re- 
cently put forward by Zhou et al |^ to interpret 
experimental data of conduction in a metal/self- 
assembled monolayer/ metal heterostructure. 

Electron tunneling rates through both metal- 
molecule junctions are obtained from a Golden rule 
calculation. The current is calculated by solving 
a master equation connecting the different charge 
states of the molecule (see Appendix A). 

To calculate the I-V characteristics, we have 
taken the apparent position of the molecular levels 
given by the T{E) spectrum of Fig. ^ in order to 
account for the shift of the level upon adsorption 
of the molecule onto the electrodes. For the sake 
of simplicity, we adjusted the experimental curve 
considering only two levels (HOMO and LUMO) 
because the charging of the molecule will most 
likely involve these two levels. The best adjust- 
ment of the experimental data is obtained assum- 
ing the HOMO to be located at the equilibrium 



Fermi level |50|. For this calculation, a charging 
energy of £"£=0.19 eV was used, which is compa- 
rable to values reported for Cgo molecules by Po- 
rath et al |Q . This is about one order of magni- 
tude smaller than the calculated charging energy 
for an isolated molecule but a strong reduction of 
the charging energy is expected for molecules ad- 
sorbed onto metallic electrodes. 
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FIG. 11. Experimental symmetric I-V curve (dia- 
monds). I-V curve (solid line) calculated within the 
framework of the sequential model. The calculation 
(see text) has been done using: Ef = Enomo, (p— 0.47 
eV and d = 0.53 nm . 

As can be seen in Fig. ^ a good qualitative de- 
scription of the experimental data can be obtained 
once the tunneling barrier parameters are adjusted 
to fit the magnitude of the current step. A prefac- 
tor is included to account for the barrier suppres- 
sion at high bias |^ (See Appendix A). The phe- 
nomenological parameters (accounting for the bar- 
rier height and width respectively) used to fit the 
data are 6 =0.47eV and d=5.3A. This approxima- 
tion is a very crude way to account for the contact 
barriers at both ends of the molecule. However, 
though the value for the barrier width is large, the 
orders of magnitude for cf) and d arc not irrealistic. 

A more microscopic approach for the trans- 
parencies of the barriers would be to use the ex- 
tended Huckel calculation to compute the overlap 
of the molecular orbitals (involved in the change of 
the total number of electrons in the molecule) with 
the Au leads and then estimate the rates from the 
broadening of these levels upon coupling. In the 
present case, the levels involved in the transport 
are the HOMOs levels. We have calculated the 
broadening of these levels using the ESQC tech- 
nique assuming a Au-S distance of 1.905 A and 
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considering each molecular orbital separately to 
avoid any interference effect between orbitals (such 
an interference effect takes place, for example, be- 
tween the HOMO and HOMO-1 orbitals in the 
T(E) spectrum in Fig. ^ and explains the non 
unity value of the corresponding peak) . We found 
a width at half maximum of 0.4eV corresponding 
to a transfer rate of about 3 10^'' Hz. This is too 
high when compared to the experimental rates of 
the order of 10^^ Hz. To get the right order of 
magnitude artificially large Au-S distances have 
to be used (typically 6 A). This shows that the 
calculation of the coupling strength based on the 
extended Hiickel model likely leads to an overes- 
timated value as already pointed out in the linear 
regime subsection above. At present, this micro- 
scopic approach unfortunately cannot be used to 
estimate the barrier suppression at high bias, be- 
cause the extended Hiickel calculation do not take 
into account the modification of the orbitals under 
bias. 



3. Discussion 

To summarize, both the coherent model at low 
bias and the sequential model at large bias give 
qualitative fits of the experimental I-V curves. 
This provides support to the hypothesis that the 
step like features in the I-V curve indeed originate 
in the discreteness of the molecular levels. How- 
ever a quantitative fit of the data has not been 
achieved. This may be explained as follows. First, 
the effect of the applied bias on the relative posi- 
tion of the resonances in the T{E) spectrum is ne- 
glected, since we consider that the bias effect is to 
shift T{E) as a whole. Although the electric fields 
involved in the experiments are always lower than 
1 X 10^ V/m, the position of the levels is expected 
to be modified |52| . Preliminary calculation using 
the DFT formalism have shown that the HOMO- 
LUMO gap of the bisthiol-bisthiophene molecule is 
indeed reduced under an electric field , and this 
effect is expected to take place also with T3. An- 
other effect of the field is to modify the coupling be- 
tween the molecule and the electrodes. This should 
show up in the T(E, V) spectrum as a modification 
of the width of the peaks and consequently in the 
I(V) as a modification of the height of the steps. 

Second, in the calculation, the molecule is as- 
sumed to be in the vacuum. In reality, given the 
experimental conditions, the molecule can be as- 
sumed to be surrounded by Argon and probably a 
few solvent molecules. We have not taken this envi- 
ronment into account in our calculations. Smooth- 
ing of the theoretical I-Vs of both models is ex- 



pected due to this polarizable surrounding. 

Additional experimental results and theoretical 
calculations are clearly needed to help deciding be- 
tween sequential and coherent tunneling or a com- 
bination of both. This includes : (i) measure- 
ments in a better controlled environment (UHV), 
and at low temperature to determine the origin 
of the broadening of the steps in the TV curve, 
(ii) calculation of the potential profile along the 
molecule in order to determine wether tunneling 
barriers are present at the end of the molecule or 
not and (iii) estimation of the tunneling times in- 
volved in the processes. Furthermore, it should be 
noted that two proposed model indeed correspond 
to extreme cases. In the coherent one, the charging 
effect is neglected, whereas in the sequential one, 
the quantum coherence is neglected. In our quali- 
tative fits, it turns out that the order of magnitude 
of the molecular level broadening (coherent model) 
and charging energy(sequential model) are compa- 
rable. Thus a proper model will very likely have 
to take into account both coherent and sequential 
transport. 

Compared to the symmetric case, the interpreta- 
tion of the asymmetric I(V) curves of Fig. ^ (curve 
a) and Fig ^ is more complicated because we can no 
longer assume that the coupling to the electrodes 
is symmetric at both ends of the molecule. First, 
these data likely correspond to the case where at 
least one end of the molecule is less coupled to 
the electrode than in the case of the symmetric 
curve since i) the zero bias conductance is smaller 
than in the symmetric case and ii) the high bias 
behaviour of the current is linear-like rather than 
exponential-like. Second, the effect of mechani- 
cally pushing on the molecule (s) induces a change 
between two different stable curves. Such an effect 
has been used recently at low bias to realise an elec- 
tromechanical amplifier based on a C60 molecule 
[|^. In the present case, two explanations of the 
observations can be put forward: either the cou- 
pling of the molecule to the electrodes is changed 
by the mechanical action, or a change in the con- 
formation of the molecule was induced. 



IV. CONCLUSION 

In this paper, we have investigated the trans- 
port properties of molecules of 2,5"-bis(acetylthio)- 
5,2',5',2"-terthienyl self-assembled in the ad- 
justable gap of a metallic break junction. We have 
observed that the I-V characteristics recorded at 
room temperature are not always symmetric with 
respect to the polarity of the applied bias and show 
two different regimes: a linear regime at low bias 
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y < 0.1 V and a highly non-linear regime with 
step-like features at higher voltage. The order 
of magnitude of the measured zero bias conduc- 
tance is comparable to the theoretical calculation 
made with the ESQC technique assuming a single 
molecule in the gap of the break junction. This 
indicates that these experiments likely involve a 
very few molecules. We have obtained qualitative 
fits of the symmetric I-V curves using two differ- 
ent models of transport, a coherent one for the low 
bias range and a sequential one for the high bias 
range, both of which explicitly involve the elec- 
tronic structure of the molecule. Finally, we have 
shown that the mechanical action of decreasing the 
gap size by 0.04 nm induces a strong modification 
of the I-V curves. 

Acknowledgments : We thank C. Bureau and 
M. Devoret for enlightening discussions, C. Blond 
and T. Risler for their contribution to the sequen- 
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APPENDIX A: THE SEQUENTIAL MODEL 

1. Description of the system 

We consider a molecule weakly coupled via tun- 
nel barriers to two electrons reservoirs (Fig.p^). 
We assume that the tranfer of an electron through 
the metal-molecule-metal junction occurs through 
subsequent tunneling events at both ends of the 
molecule 



Electrode 1 Electrode 2 

Molecule 




FIG. 12. Scheme of the metal-molecule-metal junc- 
tion, consisting of a molecule weakly coupled to elec- 
trodes via tunnel barriers (hatched) 



The reservoirs are described by a continuum of 
states ( temperature T and Fermi energy Ep, oc- 
cupation according to the Fermi-Dirac statistics 
fiE-Ep) = {l + exp{{E-EF)/kT))-'). Un- 
der bias V the electrochemical potentials of the 
electrodes are given by 

( ,, = Ep + {!- ,) eV 

[ ^1 = hp — rjeV ^ ' 

where 77 is a parameter used to describe the 
strength of coupling at both ends of the molecule. 
In the case of the symmetric coupling, we have 

7?=l/2. 

The molecule is accounted for by a set of K 
discrete molecular orbitals. The state of the 
molecule is described by its energy Ei and the 
vector (A5^, A2, . . . , A^) where A^ is the occupation 
number of the /c"* molecular orbital of energy Sh ■ 
The corresponding total number of electrons in the 
molecule is ni = 'Ylik=i'^k- energy of the 
state is 

K 

E^^Y. ^k£k + Ec {n^ - ngf (A2) 
fe=i 

where Ec is the Coulomb energy and Ug is the 
ground state number of electrons in the molecule 
plus the offset charge. 
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2. Calculation of the transfer rates 



P = ker M 



The rate of transfer of electrons through the tun- 
nel barriers are calculated by application of the 
Fermi golden rule assuming that the tunneling pro- 
cesses are elastic. 



th 



The rate F^^j" for the process connecting the i 
state and state through the transfer of one elec- 
tron from electrode k onto the molecule is given by. 



The current is finally obtained by calculating the 
net number of electrons transferred from right to 
left through either junction: 



- ji ) 



(A5) 



2tt 



2 J 



* J I 



where is the tunnel matrix element coupling 
the states i and j, and pk is the density of states 
of the electrode k, assumed constant around the 
Fermi level. 

We define in the same way connecting the 
states i and j by transferring one electron from the 
molecule onto electrode k. 



T':- ^^\T,J'p,(l~f{E,-E,)) 



Assuming that the tunneling barriers have a 
height (j) and a width d, we apply a WKB type 
approximation for the tunneling matrix element 



\Ti^j[Vb)\ a exp{-2^ d) 



where Vb is the voltage drop across a tunnel junc- 
tion. 



3. Calculation of the current 

The matrix F connecting the different states of 
the molecule is [HHl 



F = F^+ + r^- + r^+ + r^- 



(A3) 



We introduce the vector P , which i*'* compo- 
nent is the probability of having the molecule in 
state i. The evolution of the system is then de- 
scribed by 



dt 



(A4) 



definin g M by Mij=Tji i ^ j and Mii= - J2j ^ij 
equation ^Jbecomes, 



—P = MP 
dt— 



\\ ■ look for a stationary solution P assuming 
that a permanent regime rapidly sets in 



[1] A. AviRAM and M. A. Ratner, Chem. Phys. 

Lett. 29, 277 (1974). 
[2] F. L. Carter 2nd Intl Symp. Molecular Elec- 
tronic Devices (M. Dekker, New York) 149 (1982). 
[3] A. AviRAM, C. Joachim, and M. Pomerantz, 

Chem. Phys. Lett. 146, 490 (1988). 
[4] G. Lambin, M. H. Delvaux, A. Calderone, R. 

Lazzaroni, j. L. Bredas, T. C. Clarke and J. 

P. Rabe, Mol. Cryst. Liq. Cryst. 235, 75 (1993). 
[5] D. M. Cyr, B. Venkataraman, G. W. Flynn, 

A. Black and G. M. Whitesides, J. Phys. 

Chem. 100, 13747 (1996). 
[6] W. Han, E. N. Durantini, T. A. Moore, A. L. 

Moore, D. Gust, P. Rez, G. Leatherman, G. 

R. Seely, N. Tao and S. M. Lindsay, J. Phys. 

Chem. B 101, 10719 (1997). 
[7] H. Nejoh Nature 353, 640 (1991). 
[8] W. MizuTANi, M. Shigeno, K. Kajimura and 

M. Ono, Ultramicroscopy 42-44, 236 (1991). 
[9] D. PORATH and O. Millo, J. Appl. Phys. 81, 

2241 (1997). 

[10] B. Michel, G. Travaglini, H. Rohrer, C. 
Joachim and M. Amrein, Z. Phys. B 76, 99 
(1989). 

[11] X. Lu, K. W. Hipps, X. D. Wang and U. 

Marzur, j. Am. Chem. Soc. 118, 7197 (1996). 
[12] L. A. BuMM, J. J. Arnold, M. T. Cygan, T. 

D. Dunbar, T. P. Burgin, L. Jones II, D. L. 

Allara, j. M. Tour, P. S. Weiss Science 271, 

1705 (1996). 

[13] C. Dekker, S. J. Tans, B. Oberndorff, R. 
Meyer and L. C. Venema, Synth. Met. 84, 853 
(1997). 

[14] R. M. Metzger, B. Chen, U. Hopfner, M. V. 
Lakshmikantham, D. Vuillaume, T. Kawai, 
X. Wu, H. Tachibana, T. V. Hughes, H. Saku- 
RAi, J. W. Baldwin, C. Hosch, M. P. Cava, 
L. Brehmer and G. L. Ashwell, J. Am. Chem. 
Soc. 119, 10455 (1997). 

[15] S. Datta, W. Tian, S. Hong, R. Reifen- 
berger, j. I. Henderson and C. P. Kubiak, 
Phys. Rev. Lett 79, 2530 (1997). 



10 



[16] A. Dhirani, P.-H. Lin, P. Guyot-Sionnest, R. 
W. Zehner and L. R. Sita J. Chem. Phys. 106, 
5249 (1997). 

[17] C. Joachim, J. K. Gimzewski, R. R. Schlitter 

and C. Chavy, Phys. Rev. Lett. 74, 2102 (1995). 
[18] C. Joachim, and J. K. Gimzewski, Europhys. 

Lett. 30, 409 (1995). 
[19] P. Sautet and C. Joachim, Chem. Phys. Lett. 

vl85, 23 (1989). 
[20] V. MujiCA, M. Kemp, A. Roitberg and M. 

Ratner, J. Chem. Phys. 104, 7296 (1996). 
[21] V. RoussET, C. Joachim, B. Rousset and N. 

Fabre J. Phys. Ill, 5, 1985 (1995). 
[22] S. J. Tans, M. H. Devoret, H. Dai, A. Thess, 

R. E. Smalley, L. J. Geerligs and C. Dekker 

Nature. 386, 474 (1997). 
[23] T. W. Ebbesen, H. J. Lezec, H. Hiura, J. W. 

Bennett, H. F. Ghaemi and T. Thio Nature. 

382, 54 (1996). 
[24] M. A. Reed et al. Science 278, 252 (1997). 
[25] C. J. MuLLER, B.J. Vleeming, M. a. Reed, 

J.J.S. Lambda, R. Hara, L JonesII, and J.M. 

Tour Nanotechnology 7, 409 (1996). 

[26] J. M. VAN RUITENBEEK, A. AlVAREZ, L 
PiNEYRO, C. GRAHMANN, P. JOYEZ, M. H. DE- 

VORET, D. Esteve and C. Urbina, Rev. Sci. lu- 
strum. 67, 108 (1995). 

[27] Ulman, Ultrathin organic films. (Academic press, 
San Diego, 1991). 

[28] J. M. Tour et al, J. Am. Chem. Soc. 117, 9529 
(1995). 

[29] C. Lebreton Ph. D. Thesis, Paris 6 (1996). 

[30] J. G. Simmons J. Appl. Phys. 34, 1793 (1963). 

[31] R. J. P. Keijsers, O. L Shklyarevskii, J. G. 
H. Hermsen and H. van Kempen, Rev. Sci. In- 
strum. 67, 2863 (1996). 

[32] For this calculation we used the following 3s, 3p 
and 3d parameters for the S atoms (S. Alvarez, 
Extended Hiickel parameters, Barcelona University 
(1989)); Ha = -20 eV, ^ii = 2.122, Ha = -11.1 
eV, Cii = 1.827, H^i = -8 eV, = 1.5. Standard 
Hoffmann parameters were used for the C and H 
atoms. 

[33] C. Kergueris, J. P. Bourgoin and S. Palacin, 
Nanotechnology 10,8 (1999) 

[34] B. LiEDBERG, Z. Yang, I. Engquist, M. Wirde, 
U. Gelius, G. Gotz, p. Bauerle, R.M. Rum- 
mel, Ch. Ziegler and W. Gopel, J. Phys. 
Chem. B 101, 5951 (1997). 

[35] H. Sellers, A. Ulman, Y. Shnidman and J. E. 
Eilers, J. Am. Chem. Soc, 115, 9389 (1993). 

[36] The right value to use with the T3 molecule is 
not available in the litterature. A self-consistent 
approach within the DFT formalism is expected 
to provide a reliable value.C. Bureau, C. Ker- 
gueris, AND J. P. Bourgoin manuscript in 
preparation. 

[37] W. TiAN, S. Datta, S. Hong, R. Reifen- 
BERGER, J. I. Henderson and C. P. Kubiak, 



J Chem. Phys 109, 2874 (1998). 

[38] The use of a shift in energy of the molecular levels 
with regard to the Fermi level of the electrodes is a 
formal convenient way to account for a more com- 
plex reality. Indeed, the adjustment of the molec- 
ular levels to the electrodes levels should be ac- 
counted for in the density of states of the metal 
electrodes, a term we consider here as a constant 
and presently hidden in the calculation. 

[39] C. Joachim and J. F. Vinuesa, Europhys. Lett. 
33, 635 (1996). 

[40] M. Magoga and C. Joachim, Phys. Rev. B 56, 
4722 (1997). 

[41] M. BuTTiKER, Y. Imry, R. Landauer and S. 
Pinhas, Phys. Rev. B 31, 6207 (1985). 

[42] Performing the same calculation on T6, the dimer 
of T3 formed by disulhde coupling, gives G= 3 nS 
at midgap p^ . This value is also compatible with 
experimental ones. Given the experimental con- 
ditions we nevertheless think that the formation 
of T6 is unlikely. An other possibility is to con- 
sider T3 conformers. Twisting the central thio- 
phene ring from planar geometry anti-anti to pla- 
nar geometry syn-syn, G was found to vary in the 
range i. 5 — 73.5 nS at midgap. Since the calcu- 
lated energy differences between the conformers 
are varying less than 1.02 kCal/mol (about 2kT) 
at the AMI level, a conformational interconversion 
can occur even at room temperature. 

[43] E.G. Emberly and G. Kirczenow in Phys. Rev. 
B 58, 10911(1998). 

[44] Such an hypothesis was also evoked by M. Reed 
/H/ et al 

[45] S. Datta, Electronic Transport in Mesoscopic 
Systems (Cambridge University Press, Cambridge, 
1995). 

[46] C. Kergueris, J. P. Bourgoin, M. Magoga, 
C. Joachim and S. Palacin, manuscript in 
preparation. 

[47] D. Porath, Y. Levi, M. Tarabiah and O. 
MiLLO, Phys. Rev. B 56, 1 (1997). 

[48] C. Ziegler, Handbook of organic conductive 
molecules and polymers: Vol. 3. Conductive poly- 
mers: spectroscopy and physical properties. Chap. 
13 (H. S. Nalwa, 1997). 

[49] C. Zhou, M. R. Deshpande, M. A. Reed, L. 
Jones II, and J. M. Tour, Appl. Phys. Lett. 71, 
611 (1997). 

[50] A temperature of 600 K was used in the calcula- 
tion to account for a possible broadening of the 
level due to the adsorption or induced by the sur- 
rounding. 

[51] A. KOROTKOV, and Y. Nazarov, Physica B 173, 
217 (1991). 

[52] R. Lazzaroni, a. Calderone, J.L. Bredas 
and J. P. Rare, J. Chem. Phys. 107, 99 (1997). 

[53] C. Bureau, C. Kergueris and J. P. Bourgoin, 
manuscript in preparation. 

[54] C. Joachim and J. K. Gimzewski Chem. Phys. 



11 



Lett. 265, 353 (1997). 

D. Weinmann, W. Haisler and B. Kramer, 
Ann. Physik., 5, 652 (1996). 



